Voids as Alternatives to Dark Energy and the Propagation of Gamma Rays through 

the Universe. 
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We test the opacity of a void Universe to TeV energy gamma rays having obtained the extra-galactic 
background light in that Universe using a simple model and the observed constraints on the star 
formation rate history. We find that the void Universe has significantly more opacity than a ACDM 
Universe, putting it at odds with observations of BL-Lac objects. We argue that while this method 
of distinguishing between the two cosmologies contains uncertainties, it circumvents any debates 
over fine-tuning. 



Observations of type la supernovae suggest that the 
expansion of the Universe is accelerating [1| . This fact is 
usually explained by the presence of some component of 
energy density which has an equation of state w < —1/3. 
This dark energy would be very small compared to the 
density of matter at early times, becoming much larger 
relative to matter at late times, which leads to a coin- 
cidence problem: The energy density of the dark energy 
must be quite finely tuned in order to explain the shape 
of the supernova data and the density of dark energy 
and matter would be equal at some redshift z ~ 0.3, in 
other words, very close to today in the history of the Uni- 
verse when viewed in units of redshift. It can be argued 
that our living so close to a particular special time in the 
history of the Universe goes against the spirit of the cos- 
mological principle. This, combined with our complete 
lack of inspiration from particle physics as to what the 
dark energy might be, has led many researchers to seek 
alternatives. However, increasing quantities of data from 
supernovae 0, galaxy surveys |j| and the CMB Q only 
seem to make the dark energy hypothesis stronger - if we 
can trust the data, it is extremely difficult to live without 
dark energy. One notable class of examples of alterna- 
tives to dark energy are modified gravity theories, where 
we change the assumption that the Universe is described 
by general relativity on large scales Q . Another class are 
void models Q , which are the models we will focus on in 
this work. These models keep the assumption of general 
relativity but remove the assumption of isotropy and ho- 
mogeneity. The idea is that we are located close to the 
centre of a very large underdense region compared to the 
average density of the Universe. Our local underdense 
region would then act like a k = — 1 FRW Universe, the 
non-zero local fife contributing the local expansion rate 
such that if we look outwards and try to recreate the 
apparent expansion history of the Universe making the 
assumption it is homogeneous then it appears to be ac- 
celerating. The metric is therefore no longer the isotropic 
and homogeneous Robertson- Walker metric but if we as- 
sume that the void is spherical, we can use the Lemaitre- 
Tolman-Bondi (LTB) metric 



where the comoving coordinates are (r, 9, </>) and time is t. 
For matter without pressure, Einstein's equations imply 
the following constraints: 
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here a dot stands for partial derivative with respect to 
t and a prime with respect to r; p(r, t) is the energy 
density of the matter and G is Newton's constant. To 
specify the model we intend to use, we have to define 
the two arbitrary functions E(r), corresponding to the 
spatial curvature, and M(r), which is simply the mass 
integrated within a comoving radial coordinate r: 
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Both of these functions are completely defined therefore 
by the choice of an initial density profile p(r, £ltb), where 
thTB = thTB (f) refers to the beginning of the LTB evolu- 
tion and is set to a constant for simplicity. We then have 
to choose E(r) in order to match the flat FRW model at 
the boundary of the hole. This choice also ensures that 
the average density inside the hole equals the one outside, 
so that an observer situated outside the void would not 
be aware, locally, of the presence of the hole. By choos- 
ing the size of the void and the shape and magnitude 
of the Underdensity we can now fit the supernova data 
without the presence of dark energy. Each LTB model 
we fit to the data is parametrised by the initial density 
profile p(r, £ltb) defined at some early time thTB- The 
two functions E(r) and M(r) based on this profile. The 
functional form of the initial density profile is based on 
Rostov's parametrisation [3] where the density profile is 
defined as follows: 
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FIG. 1: Distribution of initial density contrast, Ai of equation 
([5]) within the void at z ~ 1100 necessary to produce voids 
which fit the supernova luminosity data. 



FIG. 2: Distribution of the size of the void, T2 of equation © 
necessary to produce voids which fit the supernova data well. 



for r < Th-, where rh is the radius of the void; and 
p(r,to) = p(to) for r > r^. Given (ri,r2), the values 
of the coefficients (Ai,A2,As) and (a,/3) are chosen so 
that p{rh,to) = p(to) and the integrated mass inside the 
hole M(r h ) = An J^ h p(r,t )r 2 dr = §7rp(i )r|. There are 
therefore three independent parameters in each fit. We 
perform fits to the supernova data by varying these initial 
conditions, the position of the solar system relative to the 
centre of the void, as well as the value of Hq in the ho- 
mogeneous part of the Universe Outside the void where 
it is constant. We do this using a Metropolis Hastings 
MCMC algorithm, the results of which for the size of the 
hole and the initial density contrast are shown in figures 
[1] and [2] respectively. It is well known that in order to fit 
the supernova data, the situation needs to be fine-tuned, 
in particular we need to be very close to the centre of the 
void for the observations to fit [§] . If we look at the size 
of the volume in which the observer has to be relative to 
the overall volume of the undcrdensity, it is typically less 
than 1% which obviously constitutes a fine tuning of our 
position in the void. Also, we know that voids do exist 
in the Universe we live in, but their size is typically only 
15-30 Mpc in radius While such voids may have a 
small effect upon the reconstruction of the dark energy 
equation of state they cannot possibly explain the 
apparent acceleration present in the supernova data. A 
quick look at figure [2] shows us that any void required to 
explain the supernova data would have to be two orders 
of magnitude larger. The power spectrum of primordial 
perturbations must have non-trivial features in it in or- 
der for such a large void on a large scale to exist with a 
non- vanishing probability 

Taking a void model seriously would therefore require 
dropping many of the assumptions that we normally 
make when we do cosmology. However, so long as the 



void model can fit the data and we only have philosoph- 
ical objections to it, it is not really clear whether we 
should disregard it outright even though it has more free 
parameters. Unlike any of the alternatives, void models 
only require general relativity and matter, so it can be 
argued that they should be given more room before we 
neglect them outright. In this work, we aim to show that 
there is another tool which can be employed to distin- 
guish between the two scenarios, namely the extragalac- 
tic background light. The continual life cycle of stars 
throughout the history of the Universe creates a back- 
ground radiation in addition to the CMB. This Extra- 
galactic Background (star) Light or EBL therefore cre- 
ates a spectrum of light from the ultra-violet to the far 
end of the infrared band, up until the point at which the 
CMB starts to dominate. Because this background light 
is made of starlight, and we are surrounded by stars, it is 
quite difficult to see. Nevertheless, there are lower lim- 
its on its magnitude across the whole range of frequencies 
due to galaxy counts and positive detections from FIRAS 
in the infrared. We model the EBL following closely the 



simplified approach of Finkc et al. [12| with a particular 
numerical implementation which lends itself to different 
expansion histories. Let us outline our way of implement- 
ing the equations, the photon density of a blackbody is 
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where e = hv/m e c 2 is the dimensionless photon energy 
in units of electron mass. The dimensionless tempera- 
ture 0(m,t*) = kBT(m, t*)/m e c 2 where T(m, i*) is the 
temperature of a star of mass m and age t„ . We neglect 
the effect of mctallicity variation over time. We want to 
know how many photons of energy e are emitted per unit 



time from a star of mass m and age i* 
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where i?(m, i*) is the radius of a star of mass m and age 
and c is the speed of light. We obtain the tempera- 
ture 6 and radius R of a star of mass m and age we 
follow Finke et al. in using the results from the pa per 
by Eggleton, Fitchctt and Tout I13| as well as the correc- 
tions outlined in the Finke work [12J. For a population 
of stars of age i* but with different masses, we need a 
weighted integral over the range of masses to include the 
effects of all the different stars. We therefore define the 
specific photon generation rate (per unit solar mass) 
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where £(m) is the initial mass function (IMF) of the 
stellar population, which we assume doesn't vary over 
time. We assume a Salpcter 1955 IMF for the stars that 
are produced, namely £( TO ) oc m -2 35 across a range of 



0.1M Q < m < m r , 



100 AIq. Here wc have 



also introduced the effects of dimming, the parameter 
/(e) is the fraction of photons which do not make it out 
of galaxies from the blue end of the spectrum due to dust. 
The energy emitted per solar mass of a particular stellar 
population created at time tf orm over its lifetime up to 
the time of the observer (e.g. today) to is 
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The dust re-radiates the absorbed energy thermally and 
modelling the multiple dust grain populations is essential 
to recreate the z = spectrum, but not to calculate the 
opacity of the Universe which mainly depends upon the 
blue end of the spectrum for TeV scale photons. The 
overall spectrum of energy density at to is finally given 
by 

eu e u(e,t ) = 

e f ° PW)] {j* P art{e,to,t)+f d art (e,t ,t)} dt (10) 

Jtetart 

where p(z) is the star formation rate in units of solar 
masses per unit time per unit volume. To calculate the 
opacity of the Universe to photons we use 
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where E y is the energy of the gamma ray at t = to 
and a(y) is the angle averaged cross section for electron 
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FIG. 3: The EBL spectrum vs. the data for the two situations 
of the best fit to the supernova data in the void model (red 
dot-dashed curve) and the ACDM case (solid black curve). 
This is plotted against actual data where the EBL light has 
been detected (blue circles) and other constraints which only 
place lower limits (green triangles). References for the data 
points are set out in the caption of figure 4 of [13]. The blue 
stellar bulge of these spectra are not fitted to this data but 
arise from the fit of equation (|12[l to the star formation history 
rate. 



positron pair production. The attenuation in photons 
(i.e. the probability of them getting to z — when emit- 
ted at an earlier time t) is then given by P = exp(— r(t)). 
We adopt the following commonly used functional form 
for the star formation rate 



bz 



p* = 



1 + (z/c) d 



(12) 



and we fit the values of a, b, c, d to fit the observational 
constraints on the star formation rate p*{t) taken from 



reference (14|. We choose one of our best fit void cos- 



mologies and compare it with ACDM. For the case of the 
void cosmology, we apply the appropriate rescaling to 
the p* data outlined in 14| to make up for the difference 
between the void cosmology and the ACDM cosmology. 
Both for the void cosmology and the ACDM cosmology, 
one obtains in this way a very satisfactory z — spec- 
trum upon tuning the temperature for the re-emission of 
absorbed light by dust, as can be seen in figure[3]We have 
therefore introduced a new data set (the star formation 
rate history) and shown that it cannot on its own dis- 
criminate between the void cosmologies and ACDM if we 
test it using the 2 = spectrum (note the EBL spectrum 
and the star formation rate data are not independent of 
each other). However by moving to a new cosmology, 
we have changed the relationship between redshift and 
time, so the physical photon density at any given red- 
shift is different which is critical for the propagation of 
TeV energy gamma rays. The photons in the EBL with 
energies up to and beyond an eV can make the Universe 
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FIG. 4: Results of our code for opacity r for the two different 
situations of a void Universe and a ACDM Universe. The 
solid Black lines correspond to the energy dependent redshift 
at which the optical depth r = 1, 2, 5, 10 moving from left to 
right. The dot-dashed red lines represent the same optical 
depths for the void cosmology. The blue dashed line is the 
constraint from observations Hal. 



opaque to TeV photons due to electron-positron pair pro- 
duction. The increased density of EBL photons in void 
Universes reduces the mean free path of TeV photons and 
we should not see them coming as abundantly from ob- 
jects at very large distances. Recent observations from 
both the HESS and MAGIC Ccrenkov gamma ray detec- 
tors have shown that TeV gamma rays are arriving from 
galaxies at truly cosmological distances, (e.g. 3C 279 at 
a redshift of z = 0.536) [ljl[l6|- In figure [4] we show our 
results vs. the constraint from the HESS telescope on the 
optical depth of the Universe to gamma rays. It is clear 
that our ACDM background light does not create prob- 
lems for the propagation of photons, although it is very 
close to coming into disagreement with the data. The 
void cosmology is starting to be in conflict with the data 
between redshifts of z = 0.15 and z = 0.4 In particular, 
this puts it at odds with the derived limits obtained by 
the MAGIC Cerenkov telescope from the spectra of the 
BL-Lac objects 1ES 1101-232 and 1ES 0347-121. This 
seems to put the void scenario at odds with the data. 
It is important to point out that here we have assumed 
a particular IMF which does not vary as a function of 
redshift, neglected the evolution of mctallicity and cho- 
sen a particular uniform dust reddening. However, it is 
clear from figurc[4]that there is a significant difference be- 
tween the void scenario and the ACDM scenario which 



we expect to remain when all these other parameters are 
altered. More data with regards to the star formation 
history of the Universe, further observations of the EBL 
and observations from the Cerenkov Telescope array will 
test this conclusion and in principle allow us to test other 
alternatives to ACDM such as modified gravity and more 
general dark energy models. 
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